Draft version February 5, 2008 

Preprint typeset using I4TgX style emulateapj v. 1 1/26/04 



X-RAY SPECTRAL VARIABILITY DURING AN OUTBURST IN VI 1 18 ORI 

M. Audard', M. Gudel", S. L. Skinner^, K. R. Briggs^, F. M. Walter'*, G. Stringfellow\ R. T. Hamilton', and 

E. F. Guinan' 

Draft version February 5, 2008 

ABSTRACT 

We present results from a multi-wavelength campaign to monitor the 2005 outburst of the low-mass young 
star VI 118 Ori. Although our campaign covers the X-ray, optical, infrared, and radio regimes, we focus in 
this Letter on the properties of the X-ray emission in VI 1 18 Ori during the first few months after the optical 
outburst. Chandra and XMM-Newton detected VI 118 Ori at three epochs in early 2005. The X-ray flux and 
luminosity stayed similar within a factor of two, and at the same level as in a pre-outburst observation in 
2002. The hydrogen column density showed no evidence for variation from its modest pre-outburst value of 
A^H ~ 3 X 10^' cm"^. However, a spectral change occurred from a dominant hot plasma 25 MK) in 2002 and 
in January 2005 to a cooler plasma (~ 8 MK) in February 2005 and in March 2005. We hypothesize that the 
hot magnetic loops high in the corona were disrupted by the closing in of the accretion disk due to the increased 
accretion rate during the outburst, whereas the lower cooler loops were probably less affected and became the 
dominant coronal component. 

Subject headings: accretion, accretion disks — stars: circumstellar matter — stars: coronae — stars: pre-main- 
sequence — stars: individual (VI 1 18 Ori) — X-rays: stars 



1. INTRODUCTION 

The origin of X-rays remains an important topic for 
young accreting stars. Although scaled-up, solar-like mag- 
netic activity is the dominant mechanism in most mod- 
erately accreting young stars (Gahm 1980; Walter & Kuhil 
ll284*'Damiani et al.Ml995; Damiani & Micel ai.l995r et c: see 
Hertout 1989; Feigelson & Montmerle 1999: lGude]|l20 04 for 
reviews), there is also evidence that accretion may play an im- 
portant role in others (Kastner et al. 2002; Stelzer & Schmitt 
|2P04; Schmitt et al. 2005). Shocks in jets may also contribute 
to the soft X -ray compo nent observed i n some young stars 
JGiidel et alJ l2005a: Kas tner et alJl2005h . Highly accreting 
stars can help us to understand how the accretion disk inter- 
acts with the stellar magnetosphere and photosphere and to 
understand the importance of accretion in the X-ray emission 
in young stars. A handful of such young stars display pow- 
erful eruptive events with flux increases in the optical regime 
of a few magnitudes. Such outbursts are thought to originate 
during a rapid increase of the disk accretion rate over a short 
peri od of time, from values of 10"^ Mq yr"' to 10"^ Mq yr"' 
(see lHartmann & Kenvon 1996) . The origin of the change in 
accretion rate is subject to debate, and could be triggered by 
close companions or thermal disk instabilities. Two classes of 
erupting young stars have emerged: FU Orionis stars (FUors), 
which display outbursts of 4 magnitudes and more that last 
several decades, and EXors (named after the prototype EX 
Lup), which display smaller outbursts ( AV = 2-3 mag) that 
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last from a few months to a few years and that may occur 
repeatedly (e.g., Herbig 1977). 

On January 10, 2005, WiUiams et alJ (12005 1) reported the 
outburst of VI 118 Ori, a low-mass Mle young EXor star 
{d = 470 pc; = O.4IM0; = 1.29/?©; P„t = 2.23 ± 
0.04 d; Lh ni > 0.25 L^jjogreffffl = 3.562; logt[yr] = 6.28; 
iHillenbrand 1997; St assun etalJil999h . VI 118 Ori has been 
known for its outbursts in the past (e.g., 19 83-84 1988-90, 
1992-94, 1997-98; see Garcia Garcia & Pars amianll2000l for 
details). For this 2005 outburst, we have started a multi- 
wavelength monitoring campaign. In this Letter, we report 
first results from our campaign, and focus on the observa- 
tions in the X-ray regime obtained with Chandra and XMM- 
Newton. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. X-Ray 

About two weeks after the lAU Circular bv lWilhams et alJ 
(2005), Chandra observed VI 118 Ori, and XMM-Newton 
subsequently observed in February and March 2005. We also 
retrieved archival data of serendipitous, pre-outburst obser- 
vations of VI 1 18 Ori by XMM-Newton in October 2001 and 
by Chandra in September 2002. The latter was published by 
Ramirez et al. (2004) who reported the detection of VI 118 
Ori (source S029). Table[3provides a log of the observations. 

The Chandra data were reduced with the CIAO 3.2.1 soft- 
ware with CALDB 3.01. For the 2005 observation, we ex- 
tracted 12 events in the 0.1-10 keV range from a circle 
with a radius of 2" centered at the position of VI 1 18 Ori on 
the ACIS-S detector For the background, we used a nearby 
10" radius circle. The scaled background contribution was 
0.52 events. For the 2002 observation, in contrast, we used 
extraction radii of 15" and 45" for VI 1 18 Ori and the nearby 
background, respectively, because of the large off-axis angle 
(0 = 8. '85) of VI 1 18 Ori on the ACIS-I detector 

The XMM-Newton data were reduced with the SAS 6.1.0 
software. For the 2005 observations, events were extracted 
using circles with radii of 20" for VI 1 18 Ori and 60" for the 
background located in a nearby source-free region. For the 



2 



Audard et al. 



2001 observation, when only MOS2 data were available*, we 
used circles of similar sizes as above at the expected posi- 
tion of VI 1 18 Ori and for the background. Although ellipses 
would better represent the shape of the point-spread-function 
at large off-axis angles where VI 118 Ori fell on the MOS2 
detector, SAS does not accurately calculate effective areas for 
such shapes. 

2.2. Optical and Near-Infrared 

We obtained optical {V, R, and /) and near-infrared (J, H, 
and K) photometry using the SMARTS^ 1 .0 m and 1.3 m tele- 
scopes. Magnitudes are measured differentially with respect 
to field stars. We have used 2MASS magnitudes to convert 
the JHK differential magnitudes to true observed magnitudes. 
For the optical bands, we have used arbitrary reference mag- 
nitudes of Vref = 10.25, Rref = 9.5, and /i-ef = 8.75 to match the 
magnitudes obtained by the Villanova group. Details about 
the SMARTS data will be provided elsewhere together with 
accompanying spectra (Stringfellow et al. 2005, in prepara- 
tion). 

Photometric coverage of VI 118 Ori from the Villanova 
University Observatory was obtained from 2005 January 28 to 
March 16. Observations were carried out in standard Bessel 
V, R, and / filters. Dark and flat field frames were collected at 
the end of each night's observations. We assigned a standard 
error of 0.05 mag for these observations, estimated from the 
signal-to-noise ratio and seeing conditions. 

2.3. Radio 

The Very Large Array ( VLA^) observed V 1 1 1 8 Ori on 2005 
January 24 from 04:56 UT to 05:56 UT in hybrid BnA con- 
figuration, using the continuum mode at 8.435 GHz (3.56 cm) 
with a total bandwidth of 100 MHz. The primary and sec- 
ondary calibrators were 3C48 and J054 1-056, respectively. 
The total on-source time was 38 minutes for VI 1 18 Ori. 

Data were edited and calibrated using AlPS^. We obtain 
a 3cr upper limit to the flux density of VI 118 Ori of 53. 6 < 
0.075 mJy, corresponding to L3 6 < 2 x 10"' ergs Hz"' s"'. 

3. RESULTS 

Figure [0 shows light curves of V 1 1 1 8 Ori in the optical 
and near-infrared bands, and in X-rays. During each of the 
X-ray observations, the X-ray flux stayed relatively constant 
and showed no strong variability due to, e.g., flares. We 
have thus used the full duration of the observations for our 
spectral analysis. We fitted the background-subtracted spec- 
tra with l-T collisional ionization equilibrium models (addi- 
tional components were not needed) based on the Astrophys- 
ical Plasma Emission Code (APEC 1.3.1:'S mith et al.ll2001h 
as implemented in XSPEC (Arnaud 1996) with a photoelec- 
tric absorption component, A'h, generally free to vary and that 
uses cross sections from Balucinska-Church & McCammon 
([l992). Ab undances are relativ e to the solar photospheric val- 
ues (iGrevesse & Sauvailll998h . Tabled summarizes the best- 

* The EPIC pn was off and the orientation of the MOS 1 detector was such 
that Vll 18 Ori fell outside the field-of-view. 

' SMARTS, the Small and Medium Aperture Research Telescope Facility, 
is a consortium of universities and research institutions that operate the small 
telescopes at Cerro Tololo under contract with AURA. 

* The National Radio Astronomy Observatory (NRAO) is a facility of the 
National Science Foundation operated under cooperative agreement by As- 
sociated Universities Inc. 

' Astronomical Image Processing System (AIPS) is a software package 
developed by NRAO. 



fit models. X-ray fluxes and luminosities with 68% confidence 
ranges. 

The observations in 2001 and 2002 are useful to determine 
the pre-outburst X-ray properties of VI 1 18 Ori and compare 
them with those in the early phases of outburst in 2005. Un- 
fortunately, VI 118 Ori was not detected in 2001 down to 
Fx < 2.3 X 10"'"^ ergs cm"^ s"' (95% Bayesian upper limit; 
[ Kraft et alJfT9 9l |) . However, it was detected in 2002 (see 
iRamfrez et alJl20 04) and its X-ray spectrum was dominated 
by a hot plasma (25 MK) and low A^h (3 x 10^' cm"^) after 
we fixed the metallicity at the outburst value, Z = 0.17 (see 
Tab.Q. 

The Chandra January 2005 observation collected too few 
counts to fit a conventional binned spectrum. We have instead 
used a maximum-likelihood fit method for unbinned data as 
described in Gudel et al. (2005b). Briefly, the method com- 
pares the distribution of the event energies with the spectral 
energy distribution of template spectra characterized by the 
plasma T, Nh, and Lx, for a fixed metallicity (we adopted 
Z = 0.17). We determined that the plasma in VI 1 18 Ori was 
hot (> 18 MK), similar to the pre-outburst 2002 data. 

We fitted the XMM-Newton EPIC pn and MOS spectra of 
February 2005 simultaneously. Whereas the X-ray flux and 
A^H did not vary much compared to January 2005 and Septem- 
ber 2002, the X-ray spectrum of VI 1 18 Ori was instead dom- 
inated by a cooler 8 MK plasma with Z = 0.17. To establish 
whether a real change in plasma temperature and A^h occurred 
between 2002 and February 2005, we performed simultane- 
ous fits of these spectra, with each parameter (i.e., T and A^h) 
in turn (a) tied to be the same value in both observations but 
free to vary, and (b) allowed to vary independently. We kept 
Z fixed to 0.17 but left the normalizations free to vary in each 
observation. F-test results confirmed that the joint fit was sig- 
nificantly improved by allowing the temperature to be differ- 
ent (F = 21.59, for a probability of < 0.001 that the improve- 
ment is by chance). On the other hand, allowing different 
A^H for a model with two temperature components was not 
significantly better (F = 0.42, for a probability of 0.52). Fig- 
ure |2] shows the observed (absorbed) X-ray fluxed spectra in 
the hard 2002 pre-outburst state and in the February 2005 soft 
outburst state. 

The March 2005 XMM-Newton data were of lower qual- 
ity than those of February 2005, due to a much higher back- 
ground level resulting from solar activity. The pn data were 
more affected than the MOS data, so we fitted only the MOS 1 
and MOS2 spectra simultaneously. We obtained cool temper- 
atures of 6-16 MK for plasma abundances between Z = 
(best-fit abundance) and Z = 0.17 (fixed value of February 
2005). As above, we performed simultaneous fits of the 
March and February 2005 spectra. F-test results confirmed 
that the fit quality did not improve when the temperature 
(F = 0.56 for a probability of 0.46) or A^h (F = 0.85 for a 
probability of 0.36) were allowed to be different in the two 
observations. However, the same comparison of the March 
2005 data with that of September 2002 found a significant 
improvement in fit quality when we allowed the temperatures 
to be different {F = 6.64 with a probability of 0.01) although 
not the A^H (F = 1.47 with a probabiHty of 0.23). Thus, the 
March 2005 plasma temperature is significantly cooler than in 
2002, while in line with the temperature observed in February 
2005. In conclusion, at some point after the initial outburst, 
a marked change in plasma temperature occurred, which has 
remained unchanged on a time scale of at least one month. 
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[3.4 ±0.3] X 10"''* ergs cm-^ s"') 
2.9!J! ? X 10"'^ ergs cm-^ s"') are 



4. DISCUSSION 

The X-ray properties of the 2005 outburst of the EXor- 
type star V1118 Ori are remarkable in view of its temper- 
ature change, which was not immediate but occurred a few 
months after the onset of the outburst in the optical and near- 
infrared regimes. Furthermore, the X-ray flux and luminosity 
of V 111 8 Ori varied by only a factor of two during the out- 
burst and compared to the serendipitous observation in 2002. 
The fluxes in the optical and near-infrared varied more signifi- 
cantly, within factors of 2- 10 (see Fig.^. The hydrogen col- 
umn density stayed relatively low, A^h ^ (1 -5) x 10 cm"^, 
throughout the observations, giving no indication for addi- 
tional absorption during the outburst. 

The disappearance of the hot plasma component, replaced 
by a dominant cool plasma in February and March 2005, is 
intriguing. There was no evidence that the high temperatures 
observed in 2002 and January 2005 were due to flares. First, 
the X-ray light curves showed no evidence for flares. Sec- 
ond, flares in magnetically active stars are accompanied by 
increases in both flux and luminosity. However, the X-ray 
fluxes, e.g., in 2002 (Fx ^ 
and in February 2005 (Fx - ^-^-o 3 
similar. In addition, the X-ray luminosity in February 2005 is 
about twice as large as the X-ray luminosity in 2002. Thus, it 
is plausible that the dramatically increased accretion rate onto 
the star is responsible, directly or indirectly, for the change in 
the plasma temperature. 

When the accretion disk closed in near the star due to the 
increase in mass accretion rate, as suggested by models for ac- 
creting binaries (Eisner & Lamb 1977; Ghosh & Lamb 1979), 
it plausibly disrupted the magnetosphere. Using the scaling 
laws of Rosner et al. ( 1978), the height of a magnetic loop, 
h, is related to the temperature at its apex (most of the emis- 
sion measure in a loop has a similar temperature) and the av- 
erage electron density, in cm"^, as h = 8.41 x IQ^T^/tie 
in cm. If the loops have a density of is 10'° cm"^, as 
is consiste nt with the den sities observed in magnetically ac- 
tive stars (Ness et al."2004), we get loop heights of 5 x 10'° 
cm for T = 25 MK, i.e., of the order of the stellar radius 
(8.6 X 10'*' cm). In contrast, cooler loops (e.g., 8 MK) have 
ten times smaller heights, i.e., they are smaller than the stel- 
lar radius and thus compact compared to the dimensions of 
the accretion disk; thus, when the accretion disk grows in- 
ward due to an increased mass accretion rate, the hotter loops 
interact first with the accreting matter and are likely to be dis- 
rupted. A new magnetic configuration then may have led pref- 



erentially to compact cool loops, thus enhancing the emission 
measure of the previously undetected cool coronal compo- 
nent. 

The above properties contrast with those of another erupt- 
ing young star, VI 647 Ori, which displayed a strong X-ray 
flux increas e coincident with the optical/near-infrared out- 
burst (Kast ner et alJl2004 . In addition, VI 647 Ori displayed 
a high A^H (~ 4 x 10^^ cm~^) which was interpreted as an in- 
dication for a dense disk wind (Grosso et al. 2005). A cool 
component was detected together with the hot component in 
VI 647 Ori, but Grosso et al. (2005) interpreted the former 
component as accretion shocks onto the stellar photosphere. 
It is unlikely that the cool component in VI 1 18 Ori originates 
from shocks because free-fall velocities of matter falling from 
the truncation radius are too slow. Even if X-rays would be 
generated in shocks, they would likely be too absorbed to be 
detected (Drake 2005). 

The VI 118 Ori observations suggest that the enhanced 
mass accretion rate during its outburst likely had an indi- 
rect impact on the star's X-ray emission. Long-term multi- 
wavelength monitoring will provide important clues to under- 
standing the importance of accretion in the X-ray emission in 
highly accreting young stars. 
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TABLE 1 

Observation LOG and X-ray Properties of VI 1 18 Ori. 



Parameter 


Oct 2001 


Sep 2002 


Jan 2005 


Feb 2005 


Mar 2005 


Satellite 


XMM-Newton 


Chandra 


Chandra 


XMM-Newton 


XMM-Newton 


ObsID 


01 12590301 


2548 


6204 


0212480301 


0212480401 




40 


48 


5 


20 


20 


Observation date (UT) 


2001 Oct 3 


2002 Sep 6-7 


2005 Jan 26 


2005 Feb 18-19 


2005 Mar 21 




01:02-12:03 


12:57-02:54 


03:07-05:04 


22:36-04:50 


16:20-22:08 


Average Julian date - 2,450,000 


2,185.8 


2,524.3 


3,396.7 


3,420.6 


3,451.3 


Net counts 


<26.3 


169.7 


11.5 


494.1" 


252.5" 


Nh (W-' cm-2) 

T (MK) 

EM (10^3 j,j„-3) 




2 7+1.2 

-OA 

25 l+'^-3 

^•■'-0.3 


1 4+3.6 

46+°^ 
^°-28 

"•■'-0.2 


4.3+1-3 
7.7+'-^ 
4 0+^-^ 


1 4+0.8 
15 9+3-9 
^•^-0.5 


z/z© 




:= 0.17 


:=0.17 


"•^'-0.06 


:= 0.17 


fx(10""ergscm-2s-i)'' 


< 2.3 


3.1-3.7 


0.9-3.5 


2.6-3.1 


5.2-6.6 


Lx (1030 ergs s"')'' 


<1.1 


1.5-1.8 


0.3-1.3 


3.2-3.8 


2.5-3.2 



Note. — The uncerlainiics arc based on 689(- Baycsian conlLclcncc ranges, whereas the upper limits lor 2001 arc 95% B ay esianupper limits. 



^Net counts for the sum of pn, MOSl, and M0S2 in February 2005 and for the sum of MOSl and M0S2 in March 2005 
''x-ray luminosity and absorbed X-ray flux in the 0.1- 10 keV range, assuming = 470 pc. 
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2000 3000 3400 3450 3500 

Julian Day - 2,450,000 

Fig. 1 . — Long-term light curves of VI 118 Ori. Notice the broken time axis with different scales for the pre-outburst epoch (1,200 days) and the outburst epoch 
(160 days). Constant zero-point magnitudes were added to the optical SMARTS data to eye-match the VUlanova data (see text). The uncertainties are plotted but 
are small compared to the day-to-day variability observed in the optical and near-infrared. The left y-axis refers to X-ray data, whereas the right y-axis refers to 
the photometric optical and near-infrared data. 
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Fig. 2. — Fluxed spectra (uncorrected for photoelectric absorption) of the September 2002 (solid) and February 2005 (dash-dotted) observations of VI 1 18 Ori. 
The soft excess (0.6-1 keV) is visible in the XMM-Newton observation during the outburst. In contrast the hard (> 2 keV) in the pre-outburst Chandra data is 
evident. 



